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Abstract. The spherical-box approach is extended to calculate the resonance 
parameters and the real part of the wave function for single particle resonances in 
a potential containing the long-range Coulomb interaction. A model potential is taken 
to demonstrate the ability and accuracy of this approach. The calculated resonance 
parameters are compared with available results from other methods. It is shown that 
in the presence of the Coulomb interaction, the spherical-box approach works well for 
not so broad resonances. In particular, for very narrow resonances, the present method 
gives resonance parameters in a very high precision. 
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1. INTRODUCTION 


The investigation of continuum and resonant states is an important subject in quantum 
physics. For the theoretical determination of resonant parameters (the energy and 
the width), in addition to many approaches which are based on the connection of 
the unbound states (including resonances and antibound or virtual states) and the 
poles of the S matrix, several bound-state-like methods have been developed, e.g., 
the complex scaling method (CSM) [1, 2, 3, 4, 5, 6, 7, 8, 9], the complex absorbing 
potential (CAP) method [10, 11], the analytical continuation in the coupling constant 
(ACCC) approach [12, 13, 14, 15, 16, 17, 18, 19, 20], and the real stabilization method 
(RSM) [21, 22, 23, 24, 25, 26, 27, 28, 29, 30]. 

In the real stabilization method, the Schrödinger equation (or the Dirac equation 
in relativistic models) of the system in question is solved in a basis [21] or a box [24] 
of finite size, thus a bound state problem is always imposed. The RSM is based on 
the fact that the energy of a “resonant” state is “stable” against the change of the 
size of the basis or the box. Many efforts have been made in order to calculate more 
efficiently resonance parameters with the RSM [24, 25, 26, 27, 28, 29]. One of the 
effective ways to do so is the so called “spherical-box approach” [24] in which the energy 
and width of a resonance in finite-range potentials are determined from the variation of 
the discrete positive energies with the radius of the box. The spherical-box approach for 
short range potential has been vigorously proved [31]. In a recent work [30], the single 
neutron resonances in atomic nuclei were investigated by combining the RSM and the 
relativistic mean field model [32, 33]. 

In many atomic, molecular and nuclear processes, the long-range Coulomb 
interactions play important roles. Therefore it is necessary to develop the spherical- 
box approach for resonances in the presence of the long-range interaction and check its 
validity and accuracy. In the present work, we extend the spherical-box approach in 
order to include the influence of the Coulomb force. 

The paper is organized as follows. In section 2 we give the formalism of the 
spherical-box approach for resonances in the presence of the Coulomb interaction. The 
numerical details, the results for a model potential and discussions are given in section 3. 
Finally in section 4 we summarize our work. 


2. Formalism 


We deal with a central field problem. The radial Schrödinger equation reads (in atomic 
units) 


1d? I(l+1) Z 
“Sa” oa Vo(r) + - V(r) = EV,(r), (1) 


where V(r) is a finite-range potential and E = k?/2. Under the box boundary condition 
the continuum is discretized and one is left with a bound state problem. When the size 
of the box R is large enough, the energies of bound states, if there are any, do not change 
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with R. In the continuum region, there are some states stable against the change of 
the size of the box, i.e., the energy of each of such states is almost constant with the 
variation of R; such stable states correspond to resonances. It should be noted that 
the spherical-box approach is essentially different from the “variable phase approach” 
where the potential in question is amputated of its part extending beyond a distance 7 
and one needs to solve a series of scattering problems for different 7 [34]. 
In the spherical-box approach, the resonance energy is determined by the stability 
condition, 
OE 


= 0, (2) 
OR? | pon, 


the corresponding box size is labeled as R, i.e., Ey = E(R). 

The width can be evaluated from the stability behavior of the positive energy state 
against the change of the box size around R [24]. When r is large enough, the finite- 
range potential Vo(r) vanishes and W;,(r) satisfies, 

1 l 
P(r) x Tr sin(kr — > — yln(2kr) + m), (3) 
T 
where m is the phase shift. —y ln(2kr) comes from the long range Coulomb interaction 
with y = Z/k. The box boundary condition Y(R) = 0 gives 
l 
kR — = — yIn(2kR) +m = ne. (4) 
Thus the derivative of the phase shift with respect to the box size reads, 


dn _ B 7 1 dE n dy 
RT- (ym) (yR) ea O 


Around an isolated resonance, the energy Æ and the phase shift m(E) satisfy the 


following relation, 


nE) = moo E) + tan (aE) © 


Under the assumption that the phase shift from the potential scattering 7 pot( Æ) varies 


slowly with respect to the box size, i.e, Om por/OR ~ 0, one derives the formula for the 
width, 


r -/7E; 
== - = (7) 
(1- Eg) (R+ 28, Eleal”) + o& ln VEER) 


Note that the above formula goes back to equation (8) in [24] when the Coulomb 
interaction is absent, i.e., Z = 0. 


3. A model problem 


In order to check the ability and the accuracy of the spherical-box approach for 
resonances in potentials containing a long-range interaction, we solve a model problem 
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and compare our results with those predicted by other methods. We choose the following 
potential, 


Z 
V(r) = Vor?e™" + E (8) 


which has been extensively investigated for Z = 0 and/or Z = —1 [35, 24, 26, 36, 37]. 
We take Vo = 7.5 in the present work. 


3.1. Numerical procedure 


The Schrödinger equation (1) with the potential (8) is solved under the box boundary 
condition by using the shooting method [38, Chapter 18] with the forth order Runge- 
Kutta algorithm. For Z = 0, the energy calculated from the shooting method converges 
with the step size decreasing and the relative error reaches to within 1078 with a step 
size ôr = 0.001. For Z = —1, a step size ôr = 0.0001 gives the relative accuracy < 1078 
in the energy. In the following we shall present results for the potential (8) with Z = —1. 
The step size dr = 0.0001 will be used in order to get a high accuracy comparable to 
the results given in the literature. 


3.2. Resonance parameters 


R (au) 


Figure 1. The energies of s wave states calculated in a spherical-box of different sizes. 
The inset shows the first avoid crossing which is zoomed in. The red dotted curve 
represents the potential (8) around the barrier. 


We first present results for s states. In figure 1 the potential V(r) given in (8) and 
the energies of s states versus the box size are presented. By examining the figure, two 
resonances can be found with energies about 1.8 and 4.0 respectively. The former is 
very narrow as is indicated by the sharp avoid crossings. The latter, however, is above 
the barrier and very broad. 

The resonance parameters evaluated at different stable regions are given in table 1 
for the very narrow resonance. For comparison, the results obtained from the complex 
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Table 1. The energy and width of the first s wave resonance evaluated at different 
avoid crossings from (2) and (7). All quantities are in atomic units. 


R 
5.0163 
7.0547 
8.7509 

10.3918 
12.0226 
13.6542 
15.2878 
16.9243 
18.5633 
20.2045 


[37] 


Ey 


1.780 531 212 
1.780 525 482 
1.780 524 837 
1.780 524 706 
1.780 524 661 
1.780 524 606 
1.780 524 620 
1.780 524 629 
1.780 524 635 
1.780 524 634 


1.780 524 536 


| T 


6.302 9x107" 
8.764 7x107" 
9.344 1x107 
9.507 7x1075 
9.553 6x10~° 
9.566 2x107" 
9.569 6x1075 
9.570 7x107" 
9.571 1x107" 
9.571 3x107 


9.571 9x107" 
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methods are also included [36, 37]. With R increasing, the resonance energy converges 


[36] | 1.780 5 | 9.58x 10-5 


to a stable value very fast. For this narrow resonance, the spherical-box approach gives 
very precise energy and width which are comparable to the exact values given in [37]. 
One achieves a seven-significant-digit accuracy for the energy at the third avoid crossing 
with R = 8.7509 and a four-significant-digit accuracy for the width at the ninth avoid 
crossing with R = 18.5633. These high precisions are very encouraging. One can then 
safely use this spherical-box approach in studying narrow resonances with computational 
efforts much less than complex calculations. 


Table 2. The energy and width of the second s wave resonance evaluated at different 
avoid crossings from (2) and (7). All quantities are in atomic units. 


R By |T 


3.4784 
4.9943 
6.2534 
7.4361 
8.5939 
9.7483 
10.9117 
12.0986 
13.3808 


4.101 765 300 
4.053 255 228 
4.021 498 436 
3.995 348 518 
3.970 907 320 
3.945 799 270 
3.917 591 013 
3.881 431 163 
3.807 212 130 


0.582 256 450 
0.804 538 485 
0.939 289 481 
1.024 872 144 
1.083 354 259 
1.101 418 144 
1.144 053 591 
1.198 720 256 
1.329 057 297 


[37] | 4.101 494 946 | 1.157 254 428 


The second s wave resonance is quite broad. 


In table 2 we list the resonance 


parameters obtained for this resonant state at different stable points R. When 
R > 13.3808, one can not find a stable region satisfying the condition 0?7E/OR? = 
0 in the present numerical accuracy. At the last observed avoid crossing, the resonance 
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energy deviates from the exact value by 8% and the width by 15% [37]. This may 
be used as a guide to estimate the accuracy of the spherical-box approach for broad 
resonances. 

A third s wave resonance was predicted to be lying at 4.66 by the complex 
method [37]. It is very broad with a width 5.34. There is no hint for this resonance 
from figure 1. This implies that such a broad resonance is beyond the ability of the 
spherical-box approach. We note that for very broad resonances, the assumption that 
the phase shift from the potential scattering m,por(Z) in (6) varies slowly with respect 
to the box size may not hold. This could be one of the reasons why the present method 
does not work well for broad resonances. 


6 


E (au) 


R (au) 


Figure 2. The energies of p wave states calculated in a spherical-box of different 
sizes. The red dotted curve represents the potential (8) together with the centrifugal 
potential. 


The energies of p states are plotted as a function of the box size R in figure 2. The 
effective potential Vog(r) = V(r) + I(l + 1)/2r? with V(r) given in (8) is also shown 
in the same figure. Only one resonance is found which is close to the barrier and the 
energy is around 3.8. This state is broader than the first but narrower than the second 
s wave resonances. The parameters for this resonance are given in table 3. Both the 
energy and the width converge well with R increasing. At the sixth stable point with 
R = 9.6945, the relative deviation of the energy from the converged value (obtained 
from the last avoid crossing at around 19.8229) is less than 0.1% and that of the width 
less than 1%. 

Figure 3 shows the E ~ R plot for the d states. There is almost no pocket in 
the effective potential. But one can still find slightly stable regions in the energy range 
4.5~5.0 in the first several E ~ R curves. Indeed we could find stable points fulfilling 
the condition 0?E/OR? = 0 at the first six avoid crossings for the d state. Thus the 
resonance parameters could be calculated and they are listed in table 4. At R = 9.3164, 
the relative deviation of the resonance energy is within 3% and that of the width is 
about 8% compared with the values obtained at the previous stable point R = 8.1438. 
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Table 3. The energy and width of the p wave resonance evaluated at different avoid 
crossings from (2) and (7). All quantities are in atomic units. 


R | & |T 


3.3342 | 3.8665 | 0.1401 
4.9740 | 3.8501 | 0.2101 
6.2430 | 3.8471 | 0.2388 
7.4205 | 3.8457 | 0.2519 
8.5641 | 3.8450 | 0.2581 
9.6945 | 3.8445 | 0.2610 
10.8201 | 3.8442 | 0.2608 
11.9441 | 3.8440 | 0.2616 
13.0681 | 3.8437 | 0.2620 
14.1925 | 3.8431 | 0.2623 
15.3174 | 3.8433 | 0.2625 
16.4430 | 3.8431 | 0.2626 
17.5691 | 3.8429 | 0.2628 
18.6958 | 3.8427 | 0.2629 
19.8229 | 3.8425 | 0.2631 


R (au) 


Figure 3. The energies of d wave states calculated in a spherical-box of different 
sizes. The red dotted curve represents the potential (8) together with the centrifugal 
potential. 


Similar to the case for the second s wave resonance, one can not find further stable 
behavior in the region R > 9.3164. 
No higher partial wave resonances are found in the present study. 


3.3. Wave functions of the resonances 


Since the schrödinger equation (1) is solved within a spherical-box, we can only get 
the real part of the wave function for each resonant states. The wave functions for the 
resonances found in the present work are shown in figure 4. The behavior of the wave 
function is consistent with the width of the resonance. For the first s wave resonance 
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Table 4. The energy and width of the d wave resonance evaluated at different avoid 
crossings from (2) and (7). All quantities are in atomic units. 


R | E, |E 


3.2624 | 4.8689 | 0.7135 
4.6913 | 4.8063 | 0.9744 
5.8924 | 4.7518 | 1.1486 
7.0239 | 4.6992 | 1.2779 
8.1438 | 4.6391 | 1.3920 
9.3164 | 4.5409 | 1.5120 


PC) 


r (au) 


Figure 4. The radial wave function V(r) of the first and the second s wave resonances 
(lower panel) and of the p and d wave resonances (upper panel). Note that when 
r > 5.0, the wave functions for the two s resonances are multiplied by 10. 


which is quite narrow, the wave function is almost completely localized inside the range 
of the finite range potential Vo(r) = 7.5 r?e~". But for the broad s resonance, the radial 
wave function inside the potential barrier is much depressed and it oscillates very much 
in the asymptotic region. For the p wave resonance, one still finds a localization feature 
in its radial wave function, but not so prominent as the case of the first s resonant state. 
The wave function of the d resonance is similar to that of the second s state (note that 
the scales in the upper and lower panels are different), except that the former vanishes 
at the origin. 

The real stabilization method is based on the fact that a resonance is more or 
less localized. Therefore the energy of a resonant state is only weakly affected by the 
variation of the size of the spherical-box. This has been shown in [24, 30] as well as in 
the present work. Next we show the convergence of the real radial wave function with 
respect to the box size. The radial wave function of the p wave resonance is shown 
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R= 

R= 

R=14.1925 
R=19.8229 ——— 


fs 


0 5 10 15 20 
r(au) 


Figure 5. The radial wave function of the p resonance within a box of different size. 
Note that the wave function is multiplied by the radial coordinate r in order to see 
clearly the asymptotic behavior. 


in figure 5. In order to see clearly the asymptotic behavior, the wave function W(r) 
is multiplied by the radial coordinate r. With R increasing, more nodes appear and 
the wave function inside the potential barrier decreases only slightly. This reveals that 
the wave function of a narrow resonance is affected little by the box size R due to the 
localization property. 


4. Summary 


The spherical box-approach, which is one of the effective implementations of the real 
stabilization methods for single particle resonances, is extended to the case in which 
a long-range force such as the Coulomb interaction plays a role. The formalism is 
presented and the numerical realization is fulfilled for this approach. 

We take a model potential (8) as an example to demonstrate the ability and the 
precision of this approach. It is shown that in the presence of the Coulomb interaction, 
the spherical-box approach still works well for narrow resonances. In particular the 
present method can give resonance parameters in a quite high precision for very narrow 
resonances. The energy and width also converge reasonably fast with the box size 
increasing. Within the spherical-box, the real wave functions of these resonances are 
obtained. The localization behavior of the radial wave function is consistent with the 
width of a resonance. The convergence of the radial wave function with the increases 
of the box size is studied, which shows that the wave function of a narrow resonance is 
also “stable” against the changes of the box size. 
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